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ABSTRACT. We present a new method for obtaining sharp asymptotics of solutions of 
the defocussing nonlinear Schrodinger (NLS) equation, based on d methods and under 
essentially minimal regularity assumptions on initial data. 



1. Introduction 

The long time behavior of solutions q(x,t) of the defocusing Nonlinear Schrodinger 
(NLS) equation 

-j, f iq t + q xx -2\q\ 2 q = 0, 

\ q(x,t = 0) = qo(x) — > as \x\ — > °° 

has been studied extensively, under varying degree of smoothness assumptions on the ini- 
tial data qo [10, 2 5]|4][7). The asymptotic behavior takes the following form: as t — > <», 
one has 

(2) q{x,t) =f- 1 / 2 a(z )e" 2/(4 ' K, ' v(zo)log(8 ' ) +^(x,f), 
where 

v(z) = -^log(l-|r(z)| 2 ), |a(z)| 2 = v(z)/2, 

and 

1 r z 7 n 

arga(z) = - / log(z- *)</(log(l - + - + argr(/v(z)) - argr(z). 

Here zo — —xj (At), Y is the gamma function and the function r is the so-called reflection 
coefficient associated to the initial data qo, as described later in this section. 

Estimates on the size of the error term S ' (x,t) depend on the smoothness assumptions on 
qo. The above asymptotic form was first obtained in IflOl . The nonlinear steepest descent 
method 1 3 1 was brought to bear on this problem in J2) (see ||5'| for a pedagogic description), 
where the authors assumed the initial data possessed high orders of smoothness and decay, 
and proved that $ (x,t) satisfied 

(3) *(«,,) = *(^) 
uniformly for all 

Early in this millennium, Deift and Zhou developed some new tools for the analysis 
of Riemann-Hilbert problems, originally aimed at studying the long time behavior of per- 
turbations of the NLS equation Q. Their methods allowed them to establish long time 
asymptotics for the NLS equation with essentially minimal assumptions on the initial data 
Q. Indeed, they showed that if the initial data is in the Sobolev space 

H 1A ={feL 2 (R):xf,f'eL 2 (m)}, 
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then the (unique, weak) solution satisfies (01 with an error term d>(x,t) that satisfies, for 
any fixed < K < 1/4, 

(4) S(x,t) = &{t-^ +K )^j . 

Recently, McLaughlin and Miller |£l O have developed a method for the asymptotic 
analysis of Riemann-Hilbert problems based on the analysis of d problems, rather than the 
asymptotic analysis of singular integrals on contours. In this paper we adapt and extend 
this method to the Riemann-Hilbert problem associated to the NLS equation. The main 
point of our work is this: by using the d approach, we avoid all delicate estimates involving 
Cauchy projection operators in L p spaces (which are central to the work in Q). In place 
of such estimates, we carry out basic estimates of double integrals which involve nothing 
more than calculus. Our result is as follows. 

Theorem 1.1. For initial data qa(x) = q{x,t = 0) in the Sobolev space 
H u ={feL 2 (R): i/,/'a 2 (l)}, 

we have as t — > °°- 

q( X ,t) = f -l/2 a(zo)el .vV(4r)-,v (zo )lo g (8,) + e 

where 

v(z) = -^log(l-|Kz)| 2 ), |«(z)| 2 = v(z)/2, 

and 

arga(z) = - / log(z-s)c/(log(l - \r(s)\ 2 )) + j +ar g r(/v(z)) -argr(z). 
n J-oo 4 

The main features of this result are threefold: 

(1) the error term is an improvement (see 0; in fact our estimate on the error is sharp) 

(2) the new d method which is used to derive it affords a considerably less technical 
proof than previous results. 

(3) the method used to establish this result is readily extended to derive a more detailed 
asymptotic expansion, beyond the leading term (see the remark at the end of the 
paper). 

The solution procedure for the nonlinear Schrodinger equation may be described as 
follows: for a given function r(z) in the Sobolev space H^' 1 = {/ € L 2 (R) : zf,f' <E 
L 2 (R) , sup gR \f{z) | < 1 }, consider the following Riemann-Hilbert problem: 

Riemann-Hilbert Problem 1. Find M = M(z) — M(z;x,t) a 2 x 2 matrix, satisfying the 
following conditions: 

{M analytic on C \ E 
M+ (z) = M_ (z)V M (z) for z e E = K with V M specified below 
M = /+(9(i) asz^oo. 

The jump matrix Vm is defined on E as follows: 

/ 1 _ |„|2 p -2it9 \ 

Vm=[ r Jul ! ). 9 = 2z 2 -4zoz = 2z 2 +xz/t. 
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Next, define 



(6) 



and then set 



Afi(*,0 = lim z(M(z)-/), 



(7) q(x,t)=2i(M 1 (x,t)) l2 . 

The fact of the matter is that q(x,t) solves the nonlinear Schrodinger equation (Q]). The con- 
nection between the initial data qo(x) and the reflection coefficient r(z) is achieved through 
the spectral and inverse-spectral theory of the associated Zakharov-Shabat differential op- 
erator 

1 \d_ ( -iq 
-1 J dx { iq 



Jzf = , 



as described, for example, in 0. It is well known that if qo S H 



1.1 



more generally this spectral transform, Stf, : H lA — > , qo ^ r, is a bijection between 
these two spaces (it turns out that M is bi-Lipschitz) ifTTI . For initial data in H l > 1 , q(x,t) 
obtained from the Riemann-Hilbert problem described above is the unique weak solution 
to the NLS equation (see Q). 



then r € H^' , and 



Recent analyses of the long-time behavior of the solution of the NLS initial value 
problem (fl~|i have involved the detailed analysis of the behavior of the solution M to the 
Riemann-Hilbert problem Q] As regularity assumptions on the initial data qo are relaxed, 
the detailed analysis becomes more involved, technically. The purpose of this manuscript 
is to carry out a complete analysis of the long-time asymptotic behavior of M under the 
assumption that r € ' , as in 0, but via a d approach which replaces very technical har- 
monic analysis involving Cauchy projection operators with very straightforward estimates 
involving some explicit two-dimensional integrals. 

2. Proof of the result 

We outline the proof of the results. We begin by outlining the solutions to two model 
RHPs which we shall need in our derivation. 



1 

-1 



2.1. A model scalar RHP. Let zq € R, and a 3 = 

satisfying: 

8 is analytic and invertible for zeC\R, 
1 as z — * °°, 
U(z) (l-kWI 2 ), z<zo 
{ 8 (z) z>zo 
The unique solution to this RHP is easily checked to be the function 

1 fzo ln(l - \r(s)\ 2 



We seek a function 8 (z) 



(8) 



S(z)- 
8+(z) 



(9) 



8(z) = exp[y(z)] = exp 



2in 



-ds 



s — z 



We will need estimates on 8{z) which we take from [7| (see Proposition (2.12) and for- 
mula (2.41) of that paper). Suppose r € L°°(R) HL 2 (R) and ||r|| L °° < p < 1. Then 



(10) 



5(z)5(z) = l, 
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(11) (l-p)a < (l-p^)J < < (1-pT* < (l-p)" 2 . 
and 

(12) |^ ±1 (^)| < 1 for ±lmz>0. 
For real z, 

(13) \8+(z)8-(z)\ — 1 and, in particular, |S(z)| = l for z > z , 

(14) |5 + (z)| = |5r 1 (z)| = (l-|r(z)| 2 )i z<zo, 
and 

_L PV r° '"gC-kWI 2 ) j 

(15) A=8+8-—ei" J -°° s-z as ) where P. V. denotes the principal value. 
Also |A| = \8+8-\ = 1, and 

(16) \\8 ± -l\\ L 2 (dz) < C -^ 



L 1 



Let X ( s ) denote the characteristic function of the interval (zq — l,Zo). Then for z € 

C\ (-oo,Z ), 

(17) Y(z) = £jlog(l - \r(s)\ 2 ) -log(l - IKzo)| 2 )Z°W(^-zo + l)} 2?r .^_ z) 

+ log(l- r(zo)r / ^~ 

= j3(z,zo) + jv(z )(l-(z-zo + l)log(z-zo+l) 
+ (z - zo) log(z - zo) + log(z -zo)), 

where v(zo) = -^log(l - |r(z )| 2 ), and 

(18) P(z,z )= r{log(l-|r( 5 )| 2 )-log(l-|r(z )|VW(*-Zo + l)} 



27ri(s-z) 

On any ray = z + e~" ? 'IR+ = {z = zo +ue~^,u > 0} with O<0<7z;or-7z;<</><O 
we find 

<^Q\ llftll < flLdlg^ 

( 19 > IIPllfli(^)< 1-p . 

where the constant c is independent of and zo (see Lemma 23.3 in [fl]). Then, by standard 
Sobolev estimates, it follows that there is a constant c independent of <p and zo such that on 



(20) 



j3 (z,zo) is continuous up to z = zo 
ll/3(.,zo)||i-(^)<c||r|| fl i,o/(l-p) 

J^(z,zo)-^(zo,zo)|<^^|z-zo| 1 / 2 . 



Notice that ( f20b also provides us with estimates that are uniformly valid in C + and 
separately. 
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2.2. A model matrix RHP. We will refer to the RHP described in this section as the 
Parabolic Cylinder RHP for reason that will soon be obvious. Given ro, a complex number 

satisfying |ro| < 1, set v = -j^log ^1 — |>"o| 2 ), and then seekP(<i;), a 2 x 2 matrix valued 

function of satisfying: 

!P is analytic for % e C \ R, 
P + ^)=P.(^V P (^,f0T^Zp. 

The contour T,p consists of four rays emanating from the origin, one in each of the four 
quadrants and described parametrically as follows: 







z,- 


= re™' 4 , r 


>0}, 


Y 11 


= 




= re 3 */ 4 , 


r>o}, 


T III 
L p 


= {<? 


■ 4 


= re^l\ 


r>o}, 


V IV 
Lp 




4 


= re 7 "/ 4 , 


r >0}. 



The jump matrix Vp is defined separately on each of the four rays, as follows: 




Figure 1 . The countour E/> and the regions Q,-, i = 1 , . . . , 6. 
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The solution to this Riemann-Hilbert problem is the explicit piecewise analytic matrix 
valued function P(%), defined below (2]|5j. We first define two auxiliary matrix valued 
functions ) and (<!;), defined in C+ and C_ respectively, as follows. For % 6 C+ 

+ / e-^Divie-**'^) e"' 4 (ft,)-' [d ? (z>_ iv ( e -*/ 4 £)) - $ D-^e-""^)] 

For % e C- 

( e«l 4 D iv (e'"%) e~™l 4 (ft,)-' (D_ iv ( e »/4|)) _ f D_ iv ( e 3*/ 4 |)] \ 

I e"/ 4 (j3 l2 )-' fe(A v (e'* /4 f)) + f Av(e' W4 |)l fr 3 ™/ 4 .D_ iv (e 3ftr / 4 f ) J' 



For a fixed complex number a, the function D a (Q is a special function, defined to be 
the unique function satisfying the parabolic cylinder equation 

^(C)+Q-^+«)a,(C) = o, 

with boundary condition 

"c« e -f 2 / 4 (i+^(r 2 )), 

for£^°°, |argC|<M 

c° e - f2/4 (i + ^(r 2 )) - (2^)'/ 2 (r(-fl))- 1 e -'r«-^ 2 /4 (i + ^(r 2 )) , 

a(g j feC— ,f<«gC<^, 

C"e-f 2 / 4 (i + ^(C~ 2 )) - (2^) 1/2 (r(-a))- 1 e -^^-«-i e ? 2 /4 (i + &($- 2 )) , 

forC^-,-^<argC<-|. 
The parameters p\\ and j3i2 are defined as follows: 

r T{-iv) 
_ (2n) l l 2 e i7l l A e- KV l 2 
fel ~ n)vr(-iv) ' 

The matrix P is defined in six non-overlapping regions: 

ft, = arg(^) e (o, |) } , H 2 = || : arg(^) € (|, ^) } , 

Q3 = {^:arg(^)efe^} ! Q 4 = arg(£) € U^fj \ , 
n 5 = {t: E (f , } , ft 6 = {^: arg(^) € (£,2*) } . 

Now P(4 ) is defined in each of these regions as follows (see FigQ]): 
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For £ e £2i : P(| ) = ¥+(4 ) e^ 2<73/4 4 ~ iV03 
For 4 G H 2 : J>(§) = e''^ /4 ^ iVCT3 



1 

-r 1 



— ? 



For 4 G £2 3 : J>(£) = »P + (4) «'« <*/ 4 f ^ l-h 

(1 '° 
o H r' 2 

For 4 G i2 5 : J>(§ ) = V - (§) e^/ 4 4 - ,VCT3 

For ^ e i2 6 : J>(§ ) = MT (§) e^/ 4 4 ^ ( J 

The fact is that P defined above solves RHP 1211 . Moreover in the notation of RHP 1211 . 
we have (see 0) 

ivr r(-iv) 



(22) P? 



^r(-/v) 



(27r) I /2 c ™/4 e -irv/2 



Note also that in the setup of RHP (f2TT >. the origin is the reference point from which the 
rays emanate. However in the following sections we will need to use the rescaling % — > 

VIf(z-zo)forzoeK- 

2.3. The RHP for NLS. Let £ = R. We seek a matrix M satisfying the following condi- 
tions: 

!M analytic on C \ E 
M+ (z) = M_ (z)V M (z) for z G E = R with V M specified below 
M = l + 0(\) asz^oo. 

The jump matrix Vm defined on E as follows: 



2 = „-2itt 



l-\r\ z -re 



VM=y\Jl t l> \ J, 0=2z 2 -4z o z. 

2.3.1. Step 1: jump matrix factorization. To the left of zo we use the factorization 
( 1 W l-|r| 2 W 1 j^jre- 21 ' 6 \ 

Vm= [ t4f^ 1 A lip ' l '\ ) =UlU0U « 



To the right of zo we use the factorization 



2.3.2. 2: Extension of r and contour deformation. We would like to deform the 
contour by opening sectors, so we need extensions off R of the off-diagonal functions 
r, f, and jE^i- We choose these extension in a way that that will make precise later. 
Let 

(1) R\ be the extension of r in the sector £2i : {z : < argz < ft/4}, 

(2) R3 be the extension of ~, r ,^ in the sector ^3 : {z : 3ft/4 < argz < n}, 
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V M = U L U Ur 



Vm = W L W R 



Figure 2. The original contour M together with the jump matrix in its 
factorized forms. 

(3) R4 be the extension of in the sector Q.4: {z: n < argz < 5n/4}, 

(4) R(, be the extension of f in the sector Q.(, : {z : 7^/4 < argz < 2%). 

Next we open sectors and define A in the regions as follows where we use the extensions 




Figure 3 . The definition of A in the different sectors £1, ■, i = 1 , . . . , 6. 

Ri defined above for Ur, Ul, Wr and Wl in the sectors. The matrix A now has the diagonal 
jump Uo = (1 — \r\ 2 ) a3 on (— °°,zo] where 



(24) 



C7 3 



1 

-1 



We can remove the jump across (— °°,zo] by defining 
(25) B=A8~ a \ 

It is easily checked algebraically that B so-defined only has jumps across the four diagonal 
rays pictured below. 



2.3.3. Specify extensions to arrive at a pure d-problem. It is now time to return to the 
extensions Ri and specify them explicitly. We seek R\ satisfying 

!R\ = r on the ray {z : z = r, r > 0} 
Ri = S\ ■= ro(z-zo)- 2iv 8 2 on the ray {z : z = re in /\ r > 0} 
\dRi\ <ci\z-zo\~ l/2 + c 2 W\ in{z: z = re ie , r > 0, < 9 < 7r/4}, 
R3 satisfying 

^3 = Yq^r on the ray {z : z = -r, r > 0} 



(27) 



^3 = / 3 := ifgji (z-zo) 2iv 8- 2 on the ray {z : z = re^l\ r > 0} 
\dRi\ <ci\z~zo\- 1/2 + c 2 \r'\ in{z: z = re ie , r>0, 3n/4 < 9 < n}, 
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Figure 4. The jump matrices for B. 



R4 satisfying 

!R A = yq^r on the ray {z : z = re 5 */ 4 , r > 0} 
R 4 = U ■= -^^{z-ZoY^d 2 on the ray {z : z = -r, r > 0} 
\dR4\<ci\z-zo\- l/2 + c 2 \r'\ in{z: z = re' e 1 r > 0, % < 9 < 5k/4}, 
and R(, satisfying 

!R(, = f on the ray {z : z = r, r > 0} 
Re=f6-= ~r (z-z ) 2 ' v 8- 2 on the ray {z : z = re liK l\ r > 0} 
\dR6\<ci\z~zo\- 1/2 +c 2 \r'\ in{z: z = re ie , r > 0, 7tt/4 < 9 < 2n}, 

The reason for requiring the above conditions will be apparent soon. 

Proposition 2.1. There exist Rj,i = 1,3,4,6 satisfying conditions f|26l)-(|29l). 

Proof. For conciseness we give the proof for R\ only; the corresponding arguments for the 
other R{ follow immediately. Define 

R\(u,v) = b(avg(u + iv))r(u) + (1 —b(axg(u + iv)))f\{u + iv) 

where 

b{9) = cos(20). 

It follows that R\ satisfies the first two conditions in d26l i immediately. Note that 

d Rl = (r-f{jdb+ b -r>, 

so 

l^il < ^o)| + \r(zo) -fi\] + c 2 |r'|. 

\z-zo\ 



Note also that 

\r-r(zo)\ = 



j\'ds < /Vll^l<IMIz?((^))-Hllk( ifl) ))<^-*l 1/2 - 



Using the expansion in (fTTJ, we also obtain that 

/, = ro (z-z )- 2 ' v 8 2 



= f exp[2fv + 2j3(zo,zo)]exp[2(j3(z,zo)-j3(z ,zo))] x 
exp[2/v((z - zo) ln(z - z ) - (z - Zo + 1) Hz -Zo + !))]• 



10 



MOMAR DIENG AND K. D. T-R MCLAUGHLIN 



We choose the free parameter ro in such a way that 

(30) r exp[2iV + 2j8(z ,Zo)] = rfo). 

Thus 

r(zo)-/i = Kzo)-Kzo)(exp[2jv((z-zo)Mz-zo)-(z-zo + l)ln(z-2o + l))+2(j3(z,zo)-j3(z 0! zo))]- 
From (JT9j and ([20j we obtain, uniformly in {z : z = re ie ,r> O,O<0<7r/4}, that 



|/3(z,z )-j3(zo,zo)| = ^(Vk^oT), 



and 

Therefore 



(z-zo)Mz-zo)|<^(Vlz-zo|) 



k(zo)-/i| = Kzo){l-exp[^(v^oT)]} 
Combining these estimates yields 

(3D i^ii< ci |z-z r 1/2 - 



-c 2 r 



a 



Recall that in the definition of the parabolic cylinder RHP (f2T~b there is a free parameter 
ro. We set 



(32) 



ro = roe y ' . 



Then by construction (conditions (f26b-(f29li). the jumps for B exactly match the the jumps 
for the parabolic cylinder RHP (LTD) . Thus we can define E — BP 1 (y8t (z — Zo)) which 
will have no jumps in the plane, but is a solution of the following problem. 

(9-problem for E. Find a 2x2 matrix-valued function E satifying: 

• E is continuous in C, 

• E — ► I as z — ► °°, 



dE = EW, 



with 



W = < 



dRie 2ite d- 2 j 
-dR 3 e- 2,,e 8 2 


N 

~dR A e 2i,e 8- 2 j 
-dR 6 e- 2,,d 8 2 


o o\ 

j 



p- 1 

y 

P- 1 



in sector £2], 
in sector D.^ , 
in sector £2zt, 
m sector Q.^ 
otherwise. 



The above problem is equivalent to the following integral equation. 

(33) E=I-~ [ [ — dA(s) 

n J J s — z 



LONG-TIME ASYMPTOTICS FOR THE NLS EQUATION VIA 5 METHODS 



1 1 



2.4. Analysis of the integral equation. The integral equation ( f33l > can be written in op- 
erator notation as follows. 



(34) 
where 

(35) 



J{E) = 



[1-7] (£)=/, 

i rr ew 



dA(s). 



KJJ S — Z 

In order to invert the integral equation, we need to show that / is small in norm. The 
following proposition gives the needed estimate. 

Proposition 2.2. There is a constant c > such that for all t > 0, the following estimate 
holds. 



(36) 



< Ct 



<cr 1 / 4 



Proof. We give the details for sector £2] only as the corresponding arguments for the other 
sectors are identical with appropriate modifications. Let H £ L°°(Oi). Then 

\HdRi8- 2 e 2M ^ 



\W)\ < 



\s-z\ 

-2 



-dA(s) 



Thus OTb yields 
where 

and 

Since r Si/ 1 ' 1 QR) we have 

141 < 

< 

Moreover 



\J(H)\<C(h+I 2 ), 



\dR ]e 2,re \ 



n, s-z 



dA(s). 



lie, 



n, s-z 



■dA(s), 



h 



-1/2 „-f«v 



•s-z 



-</A(j). 



•s-z 



- dudv 



dudv <C 

\s — z\ JQ 



JQ 



s — z 



dv. 



Thus 



s — z 



\h\<C 



i?(0H) 



o 



< 



-.dv<C 



s-z 



■ du 



1/2 



< 



I Z ((H-)) 



1/2 



'-J8I 



=c/v 



=rf v. 



Using the fact ^ffie-'P 2 ™ 1 < cf^ 1 / 4 ^ 1 / 2 , we obtain 



/o y/P-V Jo Vl 



i _--»pV 



dw < 



t l / 4 J0 y/w(\- W ) - f'/ 4 ' 
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whereas 



Hence 



=dv< 



-dw < 



p ^V^p -JoV^ ~ f 1/4 Jo VI ~ ? 1/4 



„-A 2 



-dv< 



To arrive at a similar estimate for I2, we start with the following Z/ -estimate for p > 2. 



1 



v4 



s-zo 



U>((du)) 



1 



z +v IU + IV-Z0 



P/2 



d 11 



Up 



1 



dii 



\u + iv\P/ 2 
1 

{u 2 + v 2 )Pl A 



Up 



du 



Up 



1 (l+x 2 )/V4 



dx 



Up 



< cv 



i/p-1/2 



Similarly to the L 2 -estimate above, we obtain for L q where 1/q+l/p— 1. 



1 



s — z 



< Civ-pi 1 /"- 1 



L 2 ((v,oo)) 



It follows that 

\h\<c 



r e -^ v v P -ii2\ v _«a/ q -i dv+ r e - tv > v i /P -i/2, v _pn/ q -i, 

Jo Jb 



The first integral is handled similarly as above to yield the estimate 



z-^v^P-^v-^-'dvKct-^. 



To handle the second integral, we let v = /3 + w so that the integral becomes 

r e -^+w) 2 (j8+w) i/p-i/2 w i/ ? -i rfw < r e ->™\ w y/p-u2 w uc-i d , 

Jo Jo 
Using the substitution y = \ftw then yields 

g -'vy/ P -l/2| v _0|l/ ? -l rfv < cr l/4 



Combining the previous estimates we obtain 

\h\<cr'l\ 

and the result follows. 



□ 



For t sufficiently large the integral equation (|34| | may be inverted by Neumann series. It 
follows that 



(37) 



E =1 



■(,-*) 
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2.5. Asymptotics for NLS. Recall that 

M = EP(V8t(z-z o ))8°i = (l+ 1 ^ + 0(z- 2 )^ (i + — + (( V8t(z - z ))' 2 ) )(/+■■ ■) 

for z in the sectors Q.2 and £^5 . Denote by Mi the coefficient of z~ 1 in the Laurent expansion 
ofM 

Mi 

M = I+ — +..., 

z 

so that 

P? , f 81 



M{= E\-r —j= -ri n x 

Hence the contribution to the off-diagonal entries of M\ which are relevant for the asymp- 
totics (recall[7J come from Fj° and E\. From the integral equation ( f33l> satisfied by E, we 
have 

(38) Ei = ~[[EWdA, 

and one may verify, using estimates of the type employed in section l2~4l that 

i£ii <c t -y\ 

Specifically 

\Ei\< [ [ e- ,uv \r'\dudv+ ( [ e - ,uv \z - z \ l/2 dud v = 7 3 +I 4 , 

J0 Jv JO Jzo+v 

and the Cauchy-Schwarz inequality gives 

1 /2 

h<C r ( r e- ,uv du] dv< CT 3/4 . 



Similarly, the Holder inequality for l/p + l/q=l,2<p<4 yields 

h < jT v'/P - 1/2 e -^»v du y q dv <_f_j~ v 2/p-V2 e - tv l d v < cr 3/4 

The main result Theorem (II . lb then follows from 

(39) ,(*,*) =2i^ + tf(r 3 / 4 ), 



/8f 

where (^f ) 12 is given in d22b and written explicitly in terms of r(zo) using (|32l , d30b and 
(fTSt . The identity 

I r(»v) | 2 =| r(-iV) | 2 = ^(sinh(^v))" 1 
is useful in carrying out the algebraic manipulations needed. 

Remark. It is interesting to observe that the calculations we have presented actually yield 
an asymptotic expansion (of sorts) for q(x,t). For example, instead of the integral estimates 
that follow d38l , one could use ( TTTb in d38l . yielding 

(40) E l = ^JjwdA + 0(r l ), 
which in turn implies 

(41) q (x, t )=2i¥±hz + -( [fwdA] +t?(r l ). 

V8f n \JJ J 12 
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More generally, one may use the Neumann expansion for E arising from d34l ) to obtain 
an asymptotic expansion to arbitrary order. Extraction of more detail from the subsequent 
terms in the asymptotic expansion involves an analysis of explicit double integrals built out 
of the reflection coefficient and the parabolic cylinder functions. 
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